Lie symmetries and solitons in nonlinear systems with spatially inhomogeneous 
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Using Lie group theory and canonical transformations we construct explicit solutions of nonlinear 
Schrodinger equations with spatially inhomogeneous nonlinearities. We present the general theory, 
use it to show that localized nonlinearities can support bound states with an arbitrary number 
solitons and discuss other applications of interest to the field of nonlinear matter waves. 
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Introduction.- Solitons are self- localized nonlinear 
waves which are sustained by an equilibrium between dis- 
persion and nonlinearity and appear in a great variety of 
physical contexts [1]. fn particular, these nonlinear struc- 
tures have been generated recently in ultracold atomic 
bosonic gases cooled down below the Bose-Einstein tran- 
sition temperature 0, 0, 0] ■ In those systems the effective 
nonlinear interactions are a result of the elastic two-body 
collisions between the condensed atoms. 

These interactions can be controlled by the so-called 
Feschbach resonance (FR) management J5|, which has 
been used to generate bright solitons [1, 0], induce col- 
lapse 0, etc. Recently, the control in time of the con- 
densate scattering length has been the basis for many 
theoretical proposals to obtain different types of nonlin- 
ear structures such as periodic waves ^S] , shock waves ^] , 
stabilized solitons [13], etc. 

Interactions can also be made spatially dependent by 
acting on the spatial dependence of either the magnetic 
field or the laser intensity (in the case of optical control of 
FR [111) acting on the Feschbach resonances. This possi- 
bility has motivated in the last years a strong theoretical 
interest on nonlinear phenomena in Bose-Eintein conden- 
sates (BECs) with spatially inhomogeneous interactions. 
Several phenomena have been studied in quasi-one di- 
mensional scenarios such as the emission of solitons [12] 
and the dynamics of solitons when the space modulation 
of the nonlinearity is a random [l^ . linear 15], periodic 



states with an arbitrary number solitons without any ad- 
ditional external potential, something which differs dras- 
tically from the case of inhomogeneous nonlinearities. 
Our focus will be on applications to matter waves in BEG 
problems but our ideas can also be useful in the field of 
nonlinear optical systems [l9| . 

General theory. - In this paper we will consider physical 
systems ruled by the nonlinear Schrodinger equation with 
a spatially inhomogeneous nonlinearity (INLSE), i.e. 



(1) 



where V{x) is an external potential and g{x) describes 
the spatial modulation of the nonlinearity. Stationary 
solutions of the INLSE are of the form V' — 0e^*'^* where 

-^xx + V{x)c^ + g{x)^^ = Xc^, 0(±oo) = O, (2) 

A second-order differential equation A{x,u,u' ,u") — 0, 
possesses a Lie point symmetry [g^j of the form M = 
£_{x, u)d/dx 4- ri{x, u)d/du if 



m(2)A(x,m,u',u") = 



([x,u)— -\-ri{x,u) — - 
ox ou 



A{x,u,u',u") =0. (3) 



In our case, A{x, 4>, 4>x, 4'xx) is given by Eq. ^ and the 
action of the operator M*^^^ on it leads to 



[l6| . or localized function 17]. The existence and stabil- 
ity of solutions has been studied in Ref. [13] . 

In this paper we construct general classes of nonlin- 
earity modulations and external potentials for which ex- 
plicit solutions can be constructed. To do so we will 
use Lie group theory and canonical tranformations con- 
necting problems with inhomogeneous nonlinearities with 
simpler ones having an homogeneous nonlinearity. We 
will show that localized nonlinearities can support bound 
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where 

9{x) = go6-^e-2^/.rVK.)<^^, (6a) 

c{x) = ^b'{x) + K, (6b) 

= c"(x)-b{x)V'{x)-2b'{x){V{x)-X).{6c) 

for any constant K. Eqs. ([6]) allow us to construct pairs 
{V(x), g{x)} for which a Lie point symmetry exists. 

Conservation laws and canonical transformations.- It 
is known [2l[, that the invariance of the energy is asso- 
ciated to the translational invariance. The generator of 
such a transformation is of the form M = d/dx. To use 
this fact, we define the transformation 

X - fix), U = nix)c^ (7) 

where f{x) and n{x) will be determined by requiring that 
a conservation law of energy type M = d/dX exists in 
the canonical variables. Using Eqs. ([5]) and ([7]), we get 

Jo Ks) 

n{x) = 5(a;)-l/2e~■^•/'o"l/''(")''^ (8b) 

When K = the transformations preserve the Hamilto- 
nian structure, and Eq. ^ in terms ofU — b^^^^{x) and 
X = l/b{s)ds becomes 

-^+9oU' = EU, (9) 

where E = (\ - V{x)) b{x)'^ - ^b'ix)^ + ^b{x)b"{x) is a 
constant. This means that in the new variables we ob- 
tain the nonlinear Schrddinger equation (NLSE) without 
external potential and with an homogeneous nonlinearity. 
Of course not any choice of V{x) and g{x) leads to the 
existence of a Lie symmetry or an appropriate canoni- 
cal transformation (e.g. the function b{x) must be sign 
definite for U and X to be properly defined). 

Connection between the NLSE and INLSE via the 
LSE.- We can use all the known solutions of the NLSE 
([9]), e.g. solitons, plane waves and cnoidal waves to con- 
struct solutions to Eq. ([2]). Setting K = and eliminat- 
ing c{x) in Eqs. ^ we get 

g{x) = go/bixf, (10) 

and an equation relating b{x) and V{x) 

b"'{x) - 2b{x)V'{x) + Ab'{x)\ - Ab'{x)V{x) = 0. (11) 

Although we can eliminate b{x) and obtain a nonlinear 
equation for the pairs g{x) and V{x) for which there is 
a Lie symmetry, it is more convenient to work with (jlip . 
which is a linear equation. Alternatively, we can define 
p{x) — 6^/^ (a;) and get an Ermakov-Pinney equation 




FIG. 1: [Color online] (a) Examples of exponentially localized 
nonlinearities given by Eq. (|15b|) with C\ = C2 ~ 2, C3 = 
= 1/2 (blue line), and Ci = C2 = 1.25, C3 = 0,l<; = 2 
(green line), (b) Comparison of g{x) given by Eq. (|16|l for 
(;o=aj = l,Q = 0.05 (blue line) and its harmonic approxima- 
tion (jlTp (green line), (c) Example of black soliton solution of 
Eq. ((2)1 with V = Q and an inhomogeneous nonlinearity given 
by Eq. (((13) with ui = l,a = 0.3 and go — 1. (d) Inhomoge- 
neous nonlinearity given by Eq. (|16|l with uj — 1, a = 0.3 and 
go = 1 used to calculate the black soliton solution. 

whose solutions can be constructed as 

p= {aipl + 2(3ipi^2+l^iy^^ , (13) 

with a,/3,7 constant and ^j{x) being two linearly inde- 
pendent solutions of the Schrodinger equation 

- ifixx + V{x)(p ^ Xip. (14) 

This choice leads to = AW^ with A = aj-P^ and W 
being the (constant) Wronskian W = ip'-^ip2 — (pi(p2. Thus, 
given any arbitrary solution of the linear Schrddinger 
equation (|14p we can construct solutions of the nonlinear 
spatially inhomogeneous problem Eq. (0) from the known 
solutions of Eq. Thus, using the huge amount of 

knowledge on the linear Schrodinger equation we can get 
potentials V{x) for which ipi and ip2 are known and con- 
struct b{x), the canonical transformations f{x), n{x), the 
nonlinearity g{x) and the explicit solutions (f>{x). 

Systems without external potential (V{x) = 0).- As a 
first application of our ideas let us choose V(x) = 0, then 
Eq. (fTTj) becomes b"'{x) + 46' (x) A = and its solution is 

b{x) — Ci sin wa; -I- C2 cos -f C3 (A > 0), (15a) 
b{x) = Cie'^^ + C2e-^^ + C3 (A<0), (15b) 

where lo = 2^/\X\. Using Eq. pop we see that Eq. 
(|15bp corresponds to an exponentially localized nonlin- 
earity [Fig. [Ija)] and Eq. (|15ap leads to a periodic one 

g{x) — go {1 + acosLux)'^ , (16) 

[Fig. [2b, d)]. For small a this nonlinearity is approxi- 
mately harmonic [Fig. [Hb)] 



Pxx + {X-V{x))p = E/p''. (12) 



g{x) ~ (70 (1 — Sacoswx) , a ^ 1, (17) 
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We can construct our canonical transformation by using 
Eqs. ([7]) and obtain 



tan 



a ojx 
— tan — . 
a 2 



(18) 



Using any solution of Eq. © with E = juo"^ (l — a^) 
this transformation provides solutions of Eq. ^ with 
g{x) given by For example, when ijo > we can 

obtain black soliton solutions of Eq. ([2]) of the form 




X{x) 



(19) 



where X{x) is defined by Eq. ^ [Fig. [T] (c)]. We 
emphasize that this is only a simple example of the many 
posible solutions that can be constructed in such a way. 

Concerning the case given byEq. (fTSb]) . we would like 
to discuss it in more detail since we will get an interest- 
ing phenomenon from its analysis. In order to simplify 
the following formulae (without loosing any significant 
features) we restrict ourselves to a particular choice of 
the constants f2 = l,Ci=C2 = l/2 and C3 = in Eq. 
(|15b|l . thus b{x) = coshx and Eq. ([2]) with g{x) given by 
Eq. (Uni) and A = -1/4, 



3o 



cosh X 



= 



(20) 



in terms of U and X can be written as Eq. ([9]) with 
E = —1/4 being cosX{x) = — tanhx, thus < X < tt, 
and to meet the boundary conditions (j}{±oo) — one 
has to impose U{Q) = U{tt) = 0. This means that the 
original infinite domain in Eq. (|20p is mapped into a 
bounded domain for Eq. ©. It is easy to check that 
when (70 < 0, 



U{X) = 7^ 



sn(/iX, k) 
dn(/iX, k) 



(21) 



solves Eq. © provided = 1/ [4 (l - 2fc^)] and ry^ = 
k"^ (1 - k^) I [2 |go| (1 - 2A:2)] . The function U{X) sat- 
isfies C/(0) = and in order to meet ?7(7r) = 0, the 
condition /i7r = 2nK{k) where K{k) is the elliptic inte- 

gral K{k) = (l — fc^ sin^ d</3 must hold. Thus, 

to satisfy the boundary conditions, k must choosen to 
satisfy AnK (k) \J\ — 2k'^ = tt for n = 1,2,.... It can 
be shown that for every integer number n this algebraic 
equation has only a solution fc„ which means that there 
are an infinite number of solutions of Eq. ((20)) of the form 
given by Eq. ([2T|) . Moreover, each of those solutions has 
exactly n — 1 zeroes. In Fig. [2] we plot some of them 
corresponding to rt = 1, 2, 3. These solutions can be seen 
as "bound states" of several (n) solitons with alternat- 
ing phases and their existence is remarkable. When the 
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FIG. 2: [Color online] Solutions of Eq. ^ with go = 1 
corresponding to (a) n — 1, ki = 0.634493 (b) n = 2, = 
0.690727 and (c) n = 3, feg = 0.699957. 



nonlinearity is homogeneous, g{x) = go < 0, Eq. ^ has 
only one localized solution for each A, the cosh-type soli- 
ton, in other words: there are no bound states of several 
solitons. However, when g{x) is modulated and decays 
exponentially as given in Eq. (|20p we get an infinite num- 
ber of localized solutions labelled by their finite number of 
nodes. This is a novel and interesting feature of localized 
nonlinearities. 

Systems with quadratic potentials V{x) — x^ .- Any po- 
tential for which explicit solutions of Eqs. (|14p are known 
can be used to find nontrivial nonlinearities for which so- 
lutions can be constructed. Out of many possibilities we 
discuss only an example of interest for the applications 
to nonlinear matter waves in Bose-Einstein condensates 
which is y(x) oc x^ . 

Let us choose b{x) — , which leads to a quadratic 
trapping potential V{x) — x^ and a gaussian nonlin- 
earity such as the one generated by controlling the Fes- 
chbach resonances optically using a Gaussian beam (see 
e.g. (H), thus 

gix) = 50 exp {-ix^) , V{x) = x" . (22) 

Our canonical transformation is given by X(x) = 
/o" dt exp (-^2) = (7^/2) erf a;. In this case Eq. Q 
is transformed into 

- Vxx + got/' = 0. (23) 

Note that the range of X is again finite since —^/2 < 
X < and hence, we can again construct many 

localized solutions to Eq. ([2]) starting from solutions of 
([23| which satisfy the boundary conditions C/(±y^/2) = 
0. This can be done noting that for go < and any /i the 
functions 



u'-^Hx) = 



cn{fiX, fc*) 



and 



(7(2) (X) 



sn (pX, /c») 
^/2\go\ dnifiX, fc,) 



(24) 



(25) 



with fc, = 1/^2 solve Eq. ^ and that U^'^\X) and 
C/(2)(x) vanish when ^iX = {2n + 1)K (k,) and fiX = 
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001 (Consejeria de Educacion y Ciencia de la Junta de 
Comunidades de Castilla-La Mancha, Spain). 



FIG. 3: [Color online] Different solutions of Eq. ([2| for po- 
tential and nonlinearities given by Eq. (I22|l . (a) n — 1, (b) 
n = 2 and (c) n — 3. The spatial region shown corresponds 
to a; e [—6, 6]. 



2nK (k^) correspondingly. Thus we come to an infinite 
number of solutions of Eq. (j23p under zero boundary 
conditions on the new finite interval, which correspond 
to different values of /i. Finally, localized solutions of the 
NLS equation ([2]), are given by 



with 



2nK(k,) ^j?i2 sn(6>„(a;) fc^) ^ 
\/2H9o| dn(6'„(a;), fc») 



1(2;) = riK (fc,) erf a; 



1,3,, 
2,4,, 



(26) 



(27) 



It can be shown by simple asymptotic analysis that the 
last factors in Eq. ((26|) tend to zero as x ^ icxi faster 



than exp(— x /2) and that these are indeed localized solu- 
tions of our problem as it can be seen in the ones plotted 
in Fig. [HI with different numbers of zeroes (0„(a;) pos- 
sesses n — 1 zeroes). 

In conclusion, we have used Lie symmetries and canon- 
ical transformations to construct explicit solutions of the 
nonlinear Schrodinger equation with a spatially inhomo- 
geneous nonlinearity from those of the homogeneous non- 
linear Schrodinger equation. The range of nonlinearities 
and potentials for which this can be done is very wide. 
We have studied in detail the case V — with local- 
ized and periodic nonlinearities. In the former case, we 
have used our theory to construct an infinite number of 
multi-soliton bound states, something which is not posi- 
ble in the case of spatially homogeneous nonlinearities. 
Finally, we have presented an example of physical inter- 
est (harmonic trap and gaussian nonlinearity) for which 
exact solutions can be constructed. The ideas contained 
in this paper could also be applied to study time de- 
pendent problems, higher-dimensional situations, multi- 
component systems, etc. We hope that this paper will 
stimulate further research on those topics and help to 
understand the behavior of nonlinear waves in systems 
with spatially inhomogeneous nonlinearities. 
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